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In presence of unstable dimension variability numerical solutions of chaotic systems are valid only 
for short periods of observation. For this reason, analytical results for systems that exhibit this 
phenomenon are needed. Aiming to go one step further in obtaining such results, we study the 
parametric evolution of unstable dimension variability in two coupled bungalow maps. Each of 
these maps presents intervals of linearity that define Markov partitions, which are recovered for the 
coupled system in the case of synchronization. Using such partitions we find exact results for the 
onset of unstable dimension variability and for contrast measure, which quantifies the intensity of 
the phenomenon in terms of the stability of the periodic orbits embedded in the synchronization 
subspace. 



Unstable dimension variability (UDV) is a form of non- > 
hyperbolicity in which there is no continuous splitting i 
between stable and unstable subspaces along the chaotic < 
invariant set [l[ . The variability takes place when the i 
periodic orbits, embedded in the chaotic set, have a dif- 1 
ferent number of unstable directions. This is a local phe- > 
nomenon that can influence the entire phase space, and i 
create complexity in the system H-[H- Validity of trajec- 1 
tories generated by chaotic systems that exhibit UDV is i 
guaranteed for short periods [If, which decreases as the ( 
intensity of the UDV increases jy, 0] ■ , 

The intensity of the UDV can be quantified by the em- 
bedded UPOs in a nonhyperbolic attractor Q. There are 
efficient computational methods for the analysis of these 
orbits [1, 0. However, it is a time-consuming task be- 
cause the number of orbits increases with their period, 
and in many problems it is necessary to consider very 
high periods |10l - fl2| . To avoid this problem, one con- < 
structs a model so that the UDV occurs in a transversal < 
direction to a hyperbolic attractor. The dynamics in this i 
attractor is well known, and therefore, some analytical re- < 
suits can be obtained. This type of construction allows us i 
develop tools in order to shed light on the UDV [ill (3 ■ 1 
Examples of physical problems that can be handled by 
these tools are: the effect of shadowing in the kicked 
double- rotor [TH, [r| , the be ginning of the spatial activ- ■ 
ity in the three- waves model [13, LL8|> transport properties • 
of passive inertial particles incompressible flows [13 ] , and 
the chaos synchronization in coupled map lattices [20l — 
[22l | . In some cases, the study of periodic orbits embedded 
in the synchronization subspace allows the determination 
of the global behavior of coupled chaotic maps [23| . 

The lack of accurate results hinders the understand- 
ing of the UDV. Thus, the key question that this arti- 
cle will address is the analytical calculations for systems 
that present such phenomenon. In the following pages, 



we shall consider a simple spatially extended system com- 
posed by two identical bungalow maps [24j , which are 
piecewise linear, and interacts by a diffusive coupling. 
Such a system exhibits chaos synchronization and UDV 
in the transversal direction to the synchronization sub- 
space, for certain parameters intervals (25|. Besides, this 
map presents strong chaos for the entire parameter con- 
trol interval (26|. These features allows us to study the 
parameter evolution of the UDV for arbitrary periods. 

Now, we shall consider the abovementioned map, x <— > 
f a {x), given by 



fa(x) 



1 — 2a 

2 ° d 

l-2o \ 
1=2 (1 _ 



l-3a 
l-2a 
- x) + 



in which a <G (0, 1/2} 
following property [2J 



if x e h = [0, a) 
if x G h = [a, \) 

±5i ifze/ 3 = [|,l-a) 
x) if x G I4 = [1 — a, 1] 

(1) 

This map has the 



is a parameter. 
Va, the four intervals of linearity 
of the map define four Markov partitions [27j of phase 
space ui = [jla = [0, 1] (going forward, greek indexes 
range from 1 to 4). This property allows us to study 
the symbolic and, consequently, the interval dynamics of 
the system. Therefore, we hypothesize that the Markov 
partions allows an exact result for the onset of the UDV in 
the synchronization subspace of coupled bungalow maps. 

In order to do this study, we must determine all pos- 
sible itineraries. Considering the linearity of the map in 
each interval I a and the images of its ends, 
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we obtain the graph indicated in Fig. [TJ 

Consider, now, the following matrix [26l [28|: 
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FIG. 1. Possible transitions between partitions of the bunga- 113 
low map. 114 



T = 
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T] 2 
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774 774 774 



(2) 



whose eigenvalues are given by ix,2 = {Vi i 0)/2 an d 124 
t 3i 4 = 0, in which = ^77^ +4774(772 +773). 125 
It is straightforward to apprehend that matrix ©, 126 
with all r) a = 1, represents the transfer matrix T\ - as- 127 
sociated with the graph in Fig. [1] - of the map, where 128 
the element located in the line v and column r of the 129 
71— th power, [T™]„ T , represents the number of different 130 
itineraries of size n that start in the partition I v and end 131 
in the partition I T . Therefore, the topological entropy of 132 
map ([T]) is given by logarithm of the largest eigenvalue 133 
(t\) of the matrix T\ (hx — In 2) [28| . Moreover, the 134 
invariant density of the map is given by the eigenvec- 135 
tor components associated with t\. v x = l/\/T0[2 1 1 2] T 136 

(the component v[ indicates the natural measure of then? 
a-th partition). ns 

In matrix ([2]), the rj a stands for any quantity that isi39 
constant in each interval of linearity I a and multiplicativeuo 
along a trajectory. Thus, we can use the n-th power ofoi 
the matrix @ to study the dynamical properties of thci42 
map. For example, the diagonal elements of T n providci43 
the 2™ periodic sequences of size n. m 

The trace of the matrix is directly related to its eigen-145 
values by 



trT n 



Once we know the eigenvalues of T, we can determine 
the trace of T n , whatever the value of n: 



trT™ 



2 n j 

k=0 



rj le n - k [i + {-\) n - k }. 



(3) 



In Eq. ([3]) each term in the summation is related to a 
possible symbolic sequence. Thus, if 77^ = f'{x)\ x ^i a , 
then Eq. ([3]) gives the stability coefficients spectrum 



of the n-th periodic points of the map (QJ. As an 
illustration , we have associated with the intinerary 
I1I2I4I1I2IA ■ ■ ■ Iihh a point of period 3. For this case 
the coefficient of stability is the product 771772774. 

From now on we shall examine the case of two coupled 
maps. Wc shall use the following version for the coupling: 



G 



(4) 



■JJn+lJ \VnJ \fa{Vn + s(x n ~ J/n)), 

in which 5 and e can take on different (asymmetric cou- 
pling) or equal values (symmetric coupling). If either of 
them vanishes, we obtain a master-slave coupling. In any 
instance, the dynamics, for synchronization purposes, 
will depend on their sum d = 6 + e. 

The map G keeps the unit square invariant when both 
< S < 1 and < e < 1 (we will deal only with 
these intervals). This system have the property that 
the dynamics it generates leave the straight line x = y 
of the plane invariant and, consequently, the segment 
S = {{x,y) e uj 2 I < x = y < 1}. The latter is 
often called the synchronization subspace. 

Since UDV is a local phenomenon, we shall consider 
the transversal linear stability to the synchronization 
subspace. So, wc linearize the system (j4]) and diago- 
nalize it- in the basis of the Jacobian matrix - in the 
directions U|| = [1 1] T and Uj_ = [S — e] T . The quanti- 
ties associated with the directions U| and uj_, are called 
longitudinal and transversal, respectively. 

By definition, S is nonhyperbolic if there is at least 
one periodic point embedded in that subspace whose un- 
stable dimension is different from any other point in S. 
By construction, all periodic points in the set are lon- 
gitudinally unstable. Therefore, the phenomenon occurs 
in this system only in the transversal direction and, it 
is necessary that periodic points transversely stable and 
unstable coexist with each other in S. In order to study, 
in a quantitative way, the unstable dimension variability 
of the system, we must determine the unstable dimension 
of all periodic points of the map. We must also determine 
the frequency with which a typical trajectory visits the 
neighborhood of these points. As in the synchronization 
manifold the dynamics is hyperbolic and mixing, we know 
that such frequency can be obtained by the invarariant 
density given by [2!| [30| 

P( x ) = X" lim E TTT- \T> ( 5 ) 



|A||(x,p)| ; 



where B = [x, x + Ax), and the summation extends over 
all points of period p in D, whose eigenvalues associated 
with the longitudinal direction[3l| are given by Au(x,p). 
Note that expression Q give us all possible eigenvalues 
for all points fo period p = n. It is possible calculate, 
from the number of periodic points which have the 



same eigenvalue, 
follows 1321 



For this purpose, we rewrite 6 n as 



E 



2 

w 



E 



w 



r w — r w 
V2V3 V4 



3 



Replacing in we obtain 



k=0 



w /n-k 



E 

r=0 



n — 2w 



w=0 

n-2w r w-r w 
Vl r l2 r l3 ■ 



(6) 



149 Equation (|6|) gives all information required by Eq. (|5|) . 

150 Since the system is piecewise-linear, the eigenvalues ob- 

151 taincd by ri"~ 2w ri2r]3~ r vf i can be used to determine in 

152 which partition D, these periodic points are contained. 

153 Thus, taking the 77's as the trasversal eigenvalues, we de- 

154 termine the unstable dimension of partition D. On the 

155 other hand, taking the 77's as the longitudinal eigenval- 
ue ucs, and theirs coefficients, we determine the measure, 

157 i.e. the contribution of this partition to the behavior of 

158 typical trajectories in the vicinity of the synchronization 

159 manifold. This analysis allows us to quantify the unsta- 

160 ble dimension variability. 

161 First, we determine the set of parameters for which 

162 the UDV occurs. In order to calculate the beginning of 

163 this phenomenon we evaluate the coefficients of stability 

164 of each partition. Therefore, the determination of the 

165 parameters a c and d c = S + e, which are critical for the 

166 beginning and the end of unstable dimension variability, 

167 is done by calculating the possible transversal eigenvalues 192 

168 ^ ?^?7 3 77^, with 

193 

169 ?7l = -774 = ( ) (1 - d) 



>]2 



2 a 



1 -2a 



(1-d). 



171 Simply we determine which possible combinations of 

172 77" ~ 2w if 2 ~ r 774" result the largest (sup | Aj_ | ) and the low- 

173 est (inf |Aj_|) eigenvalues, in magnitude, and evaluate the 15 

174 range of existence of the UDV like follows 



sup I Aj 



1 and inf I A 1 



f. 



(7), 



176 Since we are dealing with the magnitude of the eigen- 

. 201 

177 values, we must consider only two terms. The extremes 

, 202 

178 of the spectrum of eigenvalues are then given by: |t7i| = 

179 I (i=2) (l-d)\ and \ mm \ = \^ a {l-df\ ^. From ©~ 
wo and solving for d, we have: 



1± 



d c = 



for 0<a<l/3 



1± VsEiS ^ 1/3<o<1/2 



(8) 2 



182 The dependece of d c on a is indicated by the solid lines2io 

183 in Fig. O Note that for a = 1/3 both lines intersects, 211 

184 indicating that there is no UDV in the system for that2i2 

185 value of a. This is expected for two tent maps Iineraly2i3 

186 coupled. 214 
is? The trace T n show us the eigenvalues spectrum fonis 
las a interval of time 77, as well as the number of possiblc2i6 
189 eigenvalues. We also know that the diagonal elements of2i7 



onset of.UDV 




u 1 



FIG. 2. The UDV intensity, quantified by the contrast mea- 
sure with p — 20, as a function of the parameters a (local 
dynamics) and d (coupling strength). Solid lines given by eq. 
JSJ, indicate the transition to UDV. Dashed lines denote the 
transitions between the stabilities of the fixed points da 20th 
iteration of the map. This figure is symmetric around d = 1. 



T n are related to the stability coefficients (eigenvalues) 
of the 7i-periodic points. 

Reference [|| introduces the quantity 



- Vi(p) 



V2(P) +Ml(p) 



(9) 



called contrast measure, which quantifies the intensity of 
UDV. In Eq. ©, the quantities Hi(p) read 



/'2 



(p) = E ia k ,| Q(l A ^( fc ;p)l - 1), (n) 

' |A||(fc;p)| 



in which O(-) is the Heaviside function [34|; A\\(k;p) and 
A±(k;p) are the eigenvalue associated with the longitu- 
dinal and transversal directions to the synchronization 
subspace, respectively. These eigenvalues are calculated 
on the p-periodic point labeled by k. The summation 
extend over all fixed points of the p-th iteration of map. 
Thus, for p large enough , (ixpip) gives the probability 
of visitation of a region with unstable dimension 1 or 2 
in the p-th iteration of the map. 

Now, using what was described above, we can quan- 
tify the UDV from the coefficients in Eq. ([6]) and the Eq. 
©. We must observe that the fraction of the positive 
tranversal Lyapunov exponents [35j at p-finite time is ex- 
actly given by ^(p)- This fraction is a metric dignostic 
for UDV. If we change 0(.) by \n\A±(k;p)\ in (JTOJl and 
(|TT|) . then Hi(p) + H2{p) gives (X±(p)). So, each term 
in summation gives the contribution for the transversal 
stability of S of the respective UPO. 

Figure [5] shows, in color scale, the intensity of UDV, 
quantified by the contrast measure, in the space param- 



4 



2is ctcr. Observing the figure, we notice a large region^ 

219 limited by the solid lines, in which the system is non-238 

220 hyperbolic (C p 0). There is also a large region in239 

221 which UDV is weak. For these small values of C20 themo 

222 set of periodic orbits responsible by UDV has positive24i 

223 measure, but very small. Thus, a numerical diagnos-242 

224 tic of non-hyperbolity, as the fraction of positive finite243 

225 time Lyapunov exponent, typically cannot identify such™ 

226 regions. 245 

227 In conclusion, we have seen that the synchronized sub-246 
22s space of two coupled bungalow maps presents four in-247 

229 tervals of linearity, which define four Markov partitions248 

230 of the phase space. Since UDV does not occur in the249 

231 longitudinal direction of this subspace, we are able t025t> 

232 study analitically the symbolic and interval dynamics imsi 

233 S. Pursuant to this study we found the stability cocffi-252 

234 cients of periodic points of the dynamics in the subspace253 

235 synchronization, which in turn allowed us to write an254 

236 exact expression for the contrast measure. Thus, we es-255 



tablish analytical solutions that show the onset of UDV, 
as well as the transitions between the stability of periodic 
points in parameter space. We can use this result to iden- 
tify regions in parameter space as long as the solutions 
remain valid. 

This work has only been able to touch on the a sim- 
ple dynamical system. However, the preliminary study 
reported here has highlighted the need to explore the pos- 
sibilities of finding analytical solutions to the problem of 
UDV. As this issue involves the validity of numerical so- 
lutions is important to have exact solutions for models 
that are studied. Clearly, further studies are needed to 
understand the UDV for systems with higher dimensions 
and arbitrary elements. To carry on this research we in- 
tend to study the UDV in a coupled map lattice whose 
couplings changes over time. 

This work has been made possible thanks to the partial 
financial support from the following Brazilian research 
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333 [32] Note that the term |l + (-1)" fc J in © filters only the3« 

334 terms (n — k) which are even. 342 

335 [33] The itinerary I3I4 represents the fixed point x* — 1 — a 3 " 3 

336 of the map fl}. Exactly in this point the map is non- 344 

337 differentiable and its invariant density is discontinuous 345 
33s (except for a = 1/3). However, the density on the right 346 

339 and left of the point x* are proportional to | T73 1 and 347 

340 |»j4 1, respectively. Thus, the itinerary /3J4 represents a 348 



weighted average of the dynamics, in both sides. In the 
remaining cases, each itinerary of size n is associated with 
a point of period p = n. 
[34] We define here central directions (|Aj_| = 1) as stable 
ones. 
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